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Abstract 

Techni-fermions are added as stacks of D7-D7 techni-branes within the framework of a holo- 
graphic technicolor model that has been proposed as a realization of walking technicolor. The 
stability of the embedding of these branes is determined. When a sufficiently low bulk cut-off is 
provided the fluctuations remain small. For a longer walking region, as would be required in any 
realistic model of electroweak symmetry breaking, a larger bulk cut-off is needed and in this case 
the oscillations destabilize. 
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I. INTRODUCTION 

Walking technicolor, a non-Abelian gauge theory exhibiting nearly-conformal behavior 
over a certain finite energy range, offers the potential to provide a phenomenologically 
consistent description of dynamical electroweak symmetry breaking [H [2]. This scheme 
has been argued to be consistent with the range of values of the electroweak precision 
parameters while allowing for an extended Yukawa sector with acceptable flavor changing 
neutral currents [3j H] . Characteristic of such models are distinct regions of renormalization 
group behavior. At high energy, the gauge coupling constant exhibits asymptotic freedom. 
As the distance scale increases the coupling enters an approximate scale invariance region 
of strong coupling where it slows down ("walking" region). At still larger distances, the 
gauge coupling begins to grow rapidly until it enters the confinement region where there is 
no longer any vestige of the scale symmetry. Such renormalization group walking behavior 
was first demonstrated using the perturbative 2-loop Banks-Zaks /3-function in which a large 
number of fermions provided the necessary cancellation to produce a small [5] . In addition 
to direct field theory searches for walking /3-functions, the gauge/gravity duality provides an 
alternative means to calculate the running in non-perturbative regions of coupling constant 
space. Indeed, Nunez, Papadimitriou and Piai [6] constructed a Type IIB supergravity 
background dual to a strongly coupled M = 1 SUSY gauge theory which exhibits the above 
described running behavior. This running occurs without introducing techni-flavors and as 
such is a non-Banks-Zaks renormalization group behavior. The background is obtained by 
considering a supergravity limit of a stack of Nc D5 branes wrapping a 2-cycle resulting in 
a 10 dimensional space-time with the interval 



Here C and T are integration constant parameters and the 577(2) left-invariant coordinates 
for the S 3 are given by 




P(p) coth2p 
4 



(d6 2 + sin 9 2 dip 2 + uf + w\ 




(1) 



ui\ = cos ipd6 + sin ip sin 6d(p 
uj<i = — sin ipdO + cos tp sin 6d(p 
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W3 = dip + cos 6 dip. (2) 

The function P(p) (with P'(p) = j-P{p)) obeys a BPS obtained master differential equation 
which in this supergravity limit has the leading order small Nc / C solution 

P(p) = AC [1 + 3T (sinh 4p - 4p)] 1/3 . (3) 

Using a D5 probe brane to define an appropriate 4-dimensional 't Hooft gauge coupling 
constant [6], 

9 2 N C = N c (cosh 2p+l) 1 

8tt 2 4Csinh2p [1 + 3T(sinh4p - 4p)] x /3 ' U 

whose renormalization group behavior (in p) was shown to exhibit a walking region for 

1 < P < P* ~ il n (2/3T) >> 1 with confinement for small p < 1 and asymptotic freedom 

for p in the UV, p > p*. The walking region is seen to begin for p values for which the 

coth2p prefactor in equation Q starts to become constant. This is approximately for 

p = 1. The gauge coupling remains nearly constant until 3Tsinh4p begins to grow to 

become comparable to 1. Then the coupling constant begins to decrease in the field theory 

UV region. Thus the walking region ends approximately at 3Te 4p */2 rs 1. Since the value of 

T controls the size of the walking region its value is chosen to allow there to be a significant 

energy difference between the confinement region and the UV asymptotic freedom region. 

The renormalization group running of the gauge coupling is depicted in Fig.( [I]). 

Anguelova, et al. [7J added techni-flavors to the model by embedding N F stacks of D7 
and D7 branes in a [/-shaped geometry necessary for U(Np) x U(Np) — > U(Np) symmetry 
breakdown [9J. They employed holographic techniques in [7] in order to study the vector 
and axial vector meson contributions to the electroweak precision parameters. 

The purpose of this paper is to analyze the stability of such a D7-D7 techni-brane em- 
bedding [7J E] • Reference [7J showed that the squared masses of the vector mesons were all 
positive, while in reference [8], the squared masses of the scalar and pseudoscalar mesons 
were evaluated and again claimed to all be positive. In the present paper, we critically 
reanalyze the scalar and pseudoscalar mass spectra. While we concur that the pseudoscalar 
squared masses are positive, we find that, except for cases where the UV cut-off is such as 
to allow only a very small walking region, the lowest scalar mass squared is negative thus 
indicating that the brane embedding is unstable. 
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FIG. 1: The running 't Hooft gauge coupling constant g 2 Nc/87T 2 as a function of the logarithm of 
the renormalization scale which is related to the distance p into the bulk. Each plot is for Nc = 10 
and the upper blue curve corresponds to C = 50, T = 10~ 10 while the lower red curve corresponds 
to C = 100, T = 10~ 12 . The smaller the value of C, the stronger the gauge coupling in the walking 
region. The smaller the value of T, the greater the length of the walking region. As seen in the 
plot, the walking region extends from p ~ 1, as indicated by the vertical black dashed line, to, 
for T = 10~ 10 , p ~ p* = 5.66, as indicated by the vertical blue dashed line, and, for T = 10 -12 , 
p k, p^ = 6.81 as indicated by the vertical red dashed line. The black dotted vertical line indicates 
a typical initial walking region location corresponding to the D7-D7 overlap at po = 1.8. 

II. TECHNI-BRANE EMBEDDING 

A D7 probe brane is embedded as in [HE] according to which the brane profile is obtained 
by solving the field equations obtained from the Dirac-Born-Infeld (DBI) action of the D7 
brane with induced metric (?(§). The branes are embedded according to Table (1) with the 
D7 brane having the complementary coordinates as functions of the p coordinate, 9 = 9(p) 
and ip = (p(p). 
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Table 1. 
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The associated DBI Lagrangian is [7] 



C = L ^/-detg {8) = L (f(p) + g(p) [6%p) + sin 2 6(p) V %p)}) 1/2 

with L a constant and the subscript p denoting differentiation j- while 

sinh 2 4p 



(5) 



f(p) 

g(p) 



i 



16 5 C 6 T 2 

f(p)P(p) 
2P p (p) coth2p 



P 2 (p)P» coth^2p 



1 



-r 16 2 C 2 [1 + 3T(sinh4p-4p)] 2 / 3 
2sinh4p[l + 3T(sinh4p - 4P)] 1 / 3 
16 3 C 2 T ' 



(6) 



Note that <p is a cyclic coordinate, |^ = 0, so the field equations lead to the embedding 
equations for the D7 brane 



9 = tt/2 



dC 

dtp p 



Vaip = Po) 



(r) 



\J f{p) + g(p)v 

where the integration constant is defined by the parameter po with g(po) = 
2smh4p [i+3T^(sinh4po-4po)] = ^ where in contrast to p] no additional relation among the 
parameters po, C, and T is needed. The D7 and D7 branes overlap at p = po and their 
location on this [/-shaped embedding equation is for p > p obtained from equation 0, 



dp 



9o 



( M ) 

\g{p)J 9{p) - go 



and integrating once more to find ip = <p(p) 




(8) 



(9) 



The embedding solution is depicted in Fig. ([2]) where it is seen that the azimuthal angle <p(p) 
approaches asymptotic values <p± = ±|<p(oo)| as p goes beyond the walking region p > p*, 
while it smoothly approaches its turning point at p = po- 

For values of po and p in the walking region, the embedding profile can be found explicitly 
The functions f(p) and g(p) simplify to become 



f(p) 



45^2 ' »W ^ C 2 T 

with go = g(po) ~ e 4po /4 6 C 2 T, yielding a simple integral for ip(p) 



(10) 



f(p) 



+ / dp' 

' Po 



v / e 4(p'-p ) _ 1 J 1 



dx 



y/x^ 



±Vfe 2p0 cosh-V^ ). 



(11) 



Hence, the embedding profile can be expressed as [7] 



tanh 



cpe 



-2p 



T 



±Vl - c-^p-po). 



(12) 



The asymptotic values for the embedding angles tp± can be estimated by evaluating the 
profile at the end of the walking region 



<p ± « ±<p(p*) = ±Vfe 2po [In 2 + 2(p* - p )] . 



(13) 



For p > p** > p* in the UV region where typically p** — p* = 0.4 will be taken, the metric 
also simplifies so that the functions /(p) and g(p) become 



f(p) 



_2_\2/3 



~ -3T^ c 16p/3 



4 5 C 2 



</(p) 



4 6 C 2 T 



3 16p/3 



(14) 



With p in the walking region so that go ~ 4 |^ r , the embedding profile for large p, with 
p* < p** < p, yields the approximate relation 



<p(p) = <p(p**) + 



-e3"*e 



(15) 



A refined approximation for the asymptotic angles <p± can be obtained by expressing <p(oo) 
as 

<p(oo) = <p(p**) + <Puv = <p(P**) + Jle-l^-^e 2 ^- p0 \ (16) 

V o 

Further approximating ^(p**) using the walking region profile and a small step correction 



<p(p**) = <p(p*) + A(p, 



(17) 



where equation ^ can be used to yield 
A<p = 

~ (p** P*) 



5-0 



_#(p') - go 



(p** P*) 



S'o 



g(p**)J [g(p**) - go 



e -f(p**-p,) e 2 (>*-Po)_ 



(18) 



The various contributions to the asymptotic angles for the cases depicted in Fig. ( |2J) are 
listed in Table (2) in which p** = p* +0.40 is used and the numerical integration of equation 
Q for <p(oo) is presented as the final column entry. 
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¥>(p*) 


Ay? 


v(p**) 


Puv 


<p(°o) 


Method 


io- 10 


5.66 


0.0030791 


0.0000998 


0.0031788 


0.0000935 


0.0032723 


Approximate 


io- 10 


5.66 


0.0029884 


0.0001339 


0.0031223 


0.0000888 


0.0032111 


Numerical 


io- 12 


6.81 


0.00039208 


0.00001000 


0.00040208 


0.00000938 


0.00041146 


Approximate 


io- 12 


6.81 


0.00038305 


0.00001342 


0.00039647 


0.00000891 


0.00040538 


Numerical 



Table 2. 



III. TECHNI-BRANE INSTABILITY 

In order to determine the stability of this probe brane embedding, small oscillations about 
the D7 brane profile are considered in the (x M ,p) directions [10], [11] 



8(x,p) = n/2 + A9(x,p) 

<p(x,p) = <p(p) + Ap(x,p), (19) 

where cp(p) is the above embedding solution given by equation ([9]). Substituting these along 
with their differentials, d6{x,p) = dA9 = dx^d fl A9 + dp^A9 and dp{x,p) = d<p(p) + 
dAcp(x,p) = dp-^p(p) + dx^d^Acp + dp-^Ap, into the equation (jlj), the induced D7 brane 
metric and hence DBI action are obtained. Expanding the DBI action for the D7 brane to 
quadratic order in the fluctuations A9 and Ap yields the fluctuation Lagrangian 

£ = --a v (p)d p Apd»Ap--b v (p)(—Apj --v v (p)Ap 2 

- l -a e {p)d,A9d»A9 - h e {p) (|^) " \Mp)^\ (20) 
where the Lagrangian coefficients are given by 

Mp) = mP{p) = ^v«« 



4coth2p v //(p)+^(p)^ 4 cosh 2p 

f(p)g(p) _ [g(p)-g } 3 / 2 
UXp)+9(pW 2 y/mflpj 
P{p)^f{p)+g{p)p 2 p g{p) 



aeip) 4coth2p [g( P )-go] 

h (n) - 9 ^ - 9 M h 
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FIG. 2: Embedding functions for the D7-D7 branes. The D7 branes are embedded along the 
positive tp curves and the D7 branes are embedded along the negative <p curves. The branes 
smoothly overlap at the U-shaped tip of the embedding curves at p = po = 1.8 as more closely 
depicted in the magnified inserted plots. The numerical values of the asymptotic embedding angles 
in the upper plot are ip± = ±0.00321105 and in the lower plot are cp± = ±0.000405376. The 



dashed oblique lines are the explicit embedding profiles given by equation (12) obtained from the 
simplified metric in the walking region. 



v v (p) 

Mp) 



9(P)<P 2 P 



'f(p) 



9o 



^Jf(p) + g(pW 2 P V 9(p) y/g(p)-g 



(21) 
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Separating in the x M and p variables gives the equations d 2 A<p = M 2 A<p and d 2 A9 = M 2 A9 
where d 2 is the four-dimensional D'Alembertian, while the fluctuation equations in p take 
the form 







dp \ b{p)d W~) + a ^ M '^ ~ = 0, (22) 

where ip(p) is either the Atp or AO with their corresponding a, b, M 2 and v functions. 

The D7 and D7 branes stack separately on each side of the [/-shaped embedding profile. 
In order to distinguish these separate locations a transformation of coordinates from p to z 
is introduced with the location of the D7 branes along z > branch while the D7 branes 
are along the z < branch. The stacks smoothly overlap at z = corresponding to the 
location p = p with p an additional parameter of the model. A useful transformation 
of coordinates that separately describes the branches and at the same time converts the 
fluctuation equations to the form of one-dimensional Schrodinger equations [7] is given by 



The D7 branes lie along positive z while the D7 branes lie along negative z and the overlap 
of the D7 and D7 branes occurs at the turning point of the [/-shaped embedding, that is 
at p = Po (-2 = 0). po is an arbitrary parameter of the model along with C and T. Scaling 
each fluctuation function by its respective a and b so that 

^ = W7I < 24 > 



and substituting into the Euler-Lagrange equation (22) yields the fluctuation equation of 
the Schrodinger form 

- 0" + V0 = M 2 0, (25) 

with 

V = v/a + W + W 2 (26) 

where the prime indicates differentiation with respect to z while the superpotential W is 
given in terms of the prepotential scaling function U 

U = In {ab) 1/A = ^lnafo, (27) 



with 

__. 1 fa6V 



Inverting the coordinate transformation to find p = p(z), the potential V is now expressed 
in terms of the z coordinate, V(p(z)). 

The fluctuation spectrum values of M 2 for the techni-scalar and techni-pseudoscalar 
meson modes must in general be determined numerically for various values of po, C and T. 
Insight into the spectrum can be obtained by considering values of p in the walking region 



where the metric simplifies and / and g are approximated by equation (10). From equations 



(21) the simplified ratio a/b is found to be 



a/b = 4CT- 



, ip 



;4(po-p)) ' 



(29) 



resulting in the coordinate transformation 



-Apo 



(30) 



Scaling out the a/ CTe Apo factors and introducing the new variable ( 



'CT 



±\/e 4 ( p — 1, the prepotentials for the Aip and A9 fluctuations are obtained from the 
respective products of a and b as 



U,.,-_ 



-Ina^bp 
- In a d be - 



In 



In 



c 4 



— 1 — (C 2 + 1) 



(31) 



The superpotentials follow directly by differentiating with respect to z where now the 
prime stands for differentiation with respect to ( so that 



dz VCTe^ 



c 



Lc 2 (c 2 + 1) 
i c 



(32) 



L2 (c 2 + 1) 

The superpotential for the Aip fluctuations is singular [12] , [13] . Hence the positivity of the 
energy (mass squared) cannot be proven from the SUSY (Riccati) form of the Hamiltonian. 
The singularity in the superpotential maps the energy eigenfunctions, <f) n , out of the Hilbert 
space, so that ||(W^±ip)0 n || diverges even though \ \4> n \ \ is finite. The non-vanishing explicit 
potential vq breaks the SUSY albeit with a very small coefficient for the chosen values of T 
and po- Ignoring the SUSY breaking the A8 superpotential is non-singular, however it does 
not vanish quickly enough at large ( for there to be a normalizable ground state with zero 
energy. (For the above prepotentials and their SUSY partner potentials, the zero energy 
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ground states are of the form O = e±U — ab ±l l A which, for large £, go as C^ 1 / 2 and are 
not normalizable) . Since the walking region is of finite length the potentials will be cut- 
off at large p = p\ < p*. That is there are infinite walls at p = p\. Hence the energy 
eigenfunctions must vanish at the walls and the above "would be zero energy" ground state 
wavefunction does not satisfy this boundary condition. 

In addition to the superpotential, the potential also includes supersymmetry breaking 
pieces for the A8 fluctuations. That is, while the explicit SUSY breaking potential term for 



A(p vanishes, v v = 0, it follows from equation (21) that 



^-( T ^^k^ = (^)wTTy (33) 

Differentiating the superpotentials, the potentials for each fluctuation are found to be 



v = ^w' + W 2 = ( 1 ] \W + W 2 ] 

_ ( 1 \ 1 (C 2 + 6) 
\CTe^ ) 4((2 + 1 )2 



1 ' V V {Q 



CTe 4 P° 



V e = v e /a e + ^W' 9 + W ( 



1 



dz y y \CTe^ J [(C 2 + 1 
1 \ fl (C 2 -2) , . 1 



WL + Wl 



+ Te 4 P0 . 



CTe 4 P° ) \4((2 + 1 )2 (C 2 + ly 

1 ^ Vg(C). (34) 



CTe 4 P° 



Finally multiplying equation (25) by (CTe 4po ) and expressing all in terms of the variable £, 



the Schrodinger equations for the fluctuations take the form 

-^(0 + ^(00^(0 = ^(0 

-K(0 + Ve(C)MC) = meMO, (35) 

where the prime here again indicates differentiation with respect to ( and the squared masses 
have been scaled so that 

m\ = CTMy p0 , m 2 = CTM 2 e 4p0 . (36) 

Finally, the scaled potentials take the form 

W 4( C 2 + 1)2 
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p =1.8, T=10~ 




FIG. 3: The scaled coordinate walking region potentials for Atp, on the left with minimum 
value -3/2, and for AO, Vq on the right with maximum value 1/2. 



and are plotted in Fig.Q. This agrees with the potentials given in reference [5]. For the 
small values of T necessary for a walking region of reasonable p length, the second term in 
Vg is small compared to the first term. However it nonetheless plays a pivotal role in the 
determination of the spectrum as shown in detail in the Appendix. Note that the value of 
V e at C = is V 9 (0) = 1/2 - Te ip0 « 1/2 while 1^(0) = -3/2. 

Although the Aip superpotential is singular, its potential is non-singular. The spectrum 
for this Schrodinger equation has a negative energy ground state since J d(V^(() < and 
/ <^C I C I 7 1 V<p (01 < 00 > f° r ^ 7 < 1 [ISKIS]. Indeed the A<p potential obeys these conditions 
with 

7tt 



d(V v (() = --. (38) 

Introducing a UV cut-off so that as the walking region boundary walls move in towards 
the main region of support of the A(p potential the contribution of the infinite walls to 
the ground state energy becomes more significant and the mass squared will be positive 
indicating stability of the brane's oscillation. As detailed in the Appendix, this stability 
crossing occurs at approximately the equality of the square well energy, 7r 2 /(4^), to that of 
the absolute value of the ground state energy without walls, which is numerically determined 
as 0.56. This yields a cut-off value £a = 2.10. Thus the scalar Atp fluctuations are unstable 
for walking regions of length greater than p^ — p = 0.30. 
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The numerical determination of the brane oscillation stability is obtained as follows. 
The distance of walking from the Dl-Dl overlap at p = p to p < p* is related to the 
final renormalization scale ( according to p — Po — \ l n [C 2 + !]• Using the naive ending of 
the walking region as p* the corresponding scaled coodinate is C* = \J[e A ( p *~ Pa } — 1]. For 
Po = 1.8, the T = 10 -10 case yields p* = 5.66 corresponding to = 2, 252 while the 
T = 10 -12 case yields p* = 6.81 corresponding to C* = 22,470. Both values are far beyond 
the location of the extrema of the potentials. Thus it is reasonable to cutoff the potentials 
after a modest walking length as the wavefunction of any negative energy bound state 
(m 2 < 0) will be highly localized about the origin as are the potentials. In addition the 
extended technicolor physics is expected to be present at high energy in the walking region, 
so the model, to include this new physics, must be altered. Hence a high momentum cut-off 
will be used to indicate the necessity of this new physics entering the model [7j. 

Alternatively stated, the spectrum of Dl-Dl lowest oscillation frequencies will be ob- 
tained as the cut-off is varied. Denoting the cut-off as Ca o r correspondingly p\ = 
Po + | In + 1], the Schrodinger equations for the fluctuations are solved for scalar and 
pseudoscalar meson boundary conditions at the origin. Specifically, scalar meson fluctua- 
tions are even functions of ( such that the derivative of their wavefunction vanishes at the 
origin, <p'(( = 0) = 0. For these linearized scalar oscillations the wavefunction is normalized 
to 1 at the origin as well, <ft(( = 0) = 1. The value of the lowest oscillation frequency, m 2 or 
m 2 ), is then tuned so that the wavefunction first vanishes at the cut-off, 0(Ca) = 0. If this 
is negative, the oscillation grows exponentially in time, hence the fluctuation is not stable. 
Likewise, the pseudoscalar meson fluctuations are odd functions of ( so that their wavefunc- 
tion vanishes at the origin, = 0) = 0. Furthermore, the derivative of the wavefunction 
is normalized to 1 at the origin, <ft'(( = 0) = 1. As in the scalar meson case, the value 
of the lowest oscillation frequency, m 2 or m 2 ,, is then tuned so that the wavefunction first 
vanishes (beyond the origin) at the cut-off, </>(Ca) = 0. If this is negative, the oscillation is 
growing exponentially in time, hence the fluctuation is not stable. These lowest values of the 
mass squared for the scalar and pseudoscalar A<p and A6 mesons divided by the respective 
potentials at the origin, |V^,(0)| = 3/2 and |Vg(0)| = 1/2, are plotted as a function of the 
cut-off Ca i n Fig.( [4]). As can be seen, the A<p scalar meson is stable only for cut-off values 
below Ca = 1-5 (corresponding to p values above p = 1.8 but below the cut-off value of 
Pa = 2.095). In such a case, the cut-off confines the wavefunction to be close to the origin. 
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As the cut-off potential well walls move out beyond the majority of support of the poten- 
tial V v there effect becomes irrelevant as the wavefunction is highly localized around the 
origin and the potential exhibits one negative energy bound state with asymptotic energy 
rnt ~ —0.563 (wi^,/|V^(0)| ~ —0.375). On the other hand the A<p pseudoscalar meson and 
both AO scalar and pseudoscalar meson fluctuations are stable. As the cut-off grows, the 
mass squared values are tuned to smaller and smaller positive values for a particle in an 
infinite well. As the end of the walking region is approached the validity of the simplified 
walking region metric comes into question and the cut-off must be kept less than the walking 
region limit p\ < p*. 

To summarize the case for the scalar A<p meson, when the cut-off is of order 1, the 
spectrum is dominated by the contribution arising from the infinitely high walls so that the 
oscillation frequencies are real and the brane stable. Such a situation corresponds to the 
gauge coupling exhibiting a walking behavior only over a comparatively very small range. As 
the well walls move out past the main support of the potential, the infinite well eigenenergy 
becomes smaller and smaller so that the squared mass is dominated by the potential's bound 
state energy. This is negative and approaches a constant value as the cut-off is removed (at 
least as far as the walking region cut-off pa). The D7 brane exhibits an unstable oscillation 
mode growing exponentially in time. 

As detailed in the Appendix, the potential cannot overcome the infinite well energy for 
the pseudoscalar A<p vibration even as the cut-off is removed the squared masses are positive, 
approaching zero from above. The potential for the A9 scalar and pseudoscalar meson modes 
leads to an infinite tower of bound states as given by equation (42) in the absence of a cut-off. 
The eigenvalues are small due to the size of Te 4po << 1. As the infinite well walls move in 
towards the origin, its positive contribution to the energy eigenvalue dominates that of the 
AO well. This brane oscillation mode becomes stable at approximately where the energies 
cross, given by equation (44) with A = Te 4po ~ 1.34 x 10 -7 which leads to an extremely 
large value for the cut-off, pa > 2,000. For small cut-off values, (a ~ O(10) corresponding 
to pa ~ 3, the D7 brane AO oscillation modes are stable. 

To remove the cut-off completely, or for cut-off values beyond p*, the full metric and 
corresponding functions / and g of equation ^ must be used. This requires a full numerical 
analysis. Returning to the expressions for / and g in equation (|6]), the A<p fluctuation 
equation is transformed to a Schrodinger equation form using equation (21) for a v and 
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Lowest Scalar AO Meson Mass Squared/ Ve(0) Lowest Pseudoscalar A6 Meson Mass Squared/V#(0) 




02468I0 0241SS10 



(a (a 

FIG. 4: Lowest values of the Aip squared masses divided by the potential at the origin, mB,/|V^(0)|, 
and AO squared masses divided by the potential at the origin, rrig/\ Vg(0)\, versus the scaled coordi- 
nate cut-off value £a for the scalar and pseudoscalar mesons. The D7-D7 overlap position po = 1.8 
for all plots in the panel while the A6 potential includes the vg/ag potential with T = 10 -10 



by. The change of coordinates to z in equation (23) must be numerically evaluated as 
well as inverted numerically to find p = p(z). The case of fixed parameters C = 50 and 
T = 3.33 x 10~ 10 is studied for the D7-D7 overlap location 0.2 < po < 6.2 with no cut-off. 
The potential V v = + W' is digitized as a function of z. As in the walking region case, 
the mass squared can be scaled by the value of the potential V v at the origin for each choice 
of p , the value \V v (z — 0)| = Vq. In the walking region this was the scaled mass squared 
factor 3CTM^e Ap0 /2 as plotted in the upper left panel of Fig.( |4j). The vibration modes 
were then obtained by numerically integrating the Schrodinger equation while tuning the 
squared mass parameter to yield a normalizable wavefunction for the case of no cut-off. The 
smallest squared mass eigenvalues, scaled by the minimum of the V v potential, M^/V , for 
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FIG. 5: The solid blue curve depicts the values of the smallest techni-brane oscillation frequencies 
squared divided by the potential at the origin, M^/\V (p (z = 0)|, as a function of the D7-D7 overlap 
position po f° r fixed parameters C = 50 and T = 3.33 x 1CP 10 . The dashed red curve indicates the 
squared frequencies divided by the potential at the origin, 3CTM^e 4po /2, obtained for po in the 
walking region for large cut-off (as given in the upper left panel of Fig.( |4|). The negative values 
indicate that the lowest frequency A<p oscillation mode is unstable. 



the Aip scalar meson case are shown in Fig.( [5]). The negative eigenvalue is in agreement 
with the approximate results obtained from the walking region as detailed in the Appendix. 
Regardless of the D7-D7 overlap location, the scalar A<p meson oscillation is unstable. 



IV. CONCLUSION 

The gravity dual to a M = 1 SUSY gauge theory which exhibits an approximately 
conformal ( "walking" ) region while strongly interacting was previously introduced [6] . Stacks 
of D7-D7 techni-branes corresponding to the addition of techni-fermions in the gauge theory 
were embedded in the 10 dimensional space-time as prescribed in [7]. Fluctuations of the 
embedded branes into complementary (<p and 9) space were considered [S]. The linearized 



Euler-Lagrange equations (22) were obtained for oscillations into the complementary space 



as functions of the bulk coordinate p. A coordinate transformation (23) then resulted in the 



equations of motion 018] taking the form of one-dimensional Schrodinger equations (39 ) with 
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superpotentials W = dU jdz obtained from the prepotential scaling function U = In (afe) 1//4 . 

The scalar and pseudoscalar meson mass squared spectrum was obtained by numerically 
analyzing the Schrodinger equations. With a bulk coordinate cut-off taken to be in the 
walking region only, 1 < p < p\ < P*, the fluctuation potentials were found explicitly, (cf. 
Eq. (37)), and plotted in Fig.(|3]). The A9 meson modes form an infinite tower of negative 
states in the absence of cut-off walls. As the walls are brought in from infinity these energies 
however become positive when the contribution of the square well wall's energy exceeds the 
potential's bound state energy. This occurs at extremely large values of the cut-off due to 
the smallness of T. Correspondingly the brane oscillations are stable as depicted in the two 
lower panels in Fig.( [1]). As seen in the Appendix, the mass squared eigenvalue for the A(p 
pseudoscalar meson vanishes in the SUSY (A — > 0) limit with no cut-off. When the cut-off 
is introduced all energies are positive, as depicted in the upper right panel in Fig.Q, and 
hence exhibit stable oscillations. 

On the other hand, the Atp scalar meson ground state energy arises from two sources. 
One is the negative energy contribution arising from the localized potential and the other is 
a positive term coming from the infinitely high walls from the cut-off. Unless the walls are 
very close to the position of the highly localized potential, the negative energy contribution 
of the potential dominates. Only for small values (< (a = 2.095) of the cut-off does the 
positive energy contribution of the walls takes over leading to a stable brane. As the cut-off 
increases, the positive energy contribution of the infinite walls decreases and the negative 
energy of the brane dominates the fluctuation mass squared eigenvalue, hence destabilizing 
the oscillations of the brane. As the cut-off walls move beyond the major support of the 
brane's potential, the square well contribution to the mass squared energy is negligible. The 
brane's wavefunction is highly localized in this bound state about the origin, i.e. close to po- 
The negative mass squared then attains its constant value as seen in the upper left panel of 
Fig. ( |4]) . For the no cut-off case, the full metric must be used not just its simplified form in 
the walking region to find the Schrodinger equation potentials. The numerical analysis for 
the A<p scalar meson fluctuation mass squared eigenvalue was performed as a function of 
the D7-D7 branes overlap position p . The brane oscillations are unstable. In the walking 
region the scalar meson mass squared value agreed with that calculated using the simplified 
walking region metric as depicted in Fig. ^ . 

The range of stability of the techni-brane oscillations with a cut-off in the walking region 
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is limited to 1.8 < p < 2.095. Whether a different embedding of the stacks of D7-D7 branes 
involving the 2-cycle £2 used to wrap the stacked D5 branes of the gauge theory made from 
the twisting of the S 2 coordinates 9, ip and the S 3 coordinates 9, (p and ip will lead to a 
stable techni-brane embedding remains to be investigated. 

In reference [8], the authors compared the scalar meson spectrum to that of the vector 
meson spectrum |7] and concluded the mass splitting was small and hence the brane oscil- 
lations were stable. It has been shown in this paper that the scalar meson spectrum (in 
particular that of A<p) is significantly modified compared to the vector meson case. Indeed 
as shown above there is a deep lying narrow bound state in the A(p scalar meson case with 
energy given by equation (47) that is not of negligible size. It is the presence of this state, 
whose ramifications were not properly accounted for previously [SJ, which is responsible for 
the destabilization of the brane except for the case of very small cut-offs. 
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Appendix 

The structure of the scalar and pseudoscalar mass spectra can be understood in terms 
of the approximate supersymmetry [12j [13] and conformal symmetry [T4"l [15] of equivalent 
quantum mechanical models with Hamiltonian H = — + V((). 

Model 1: The normal mode equation describing the A9 fluctuations is equivalent 
to the Schrodinger equation of a quantum mechanical model with potential 

,/(c) = -i#f^Vrr (39) 

where the potential parameter A is identified with Te 4p0 and the energy eigenvalues with the 
mass-squared values m^: Eq = m 2 B . The characteristic shape of the potential with a hump 
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at the origin flanked by two valleys is depicted in the left panel of Fig.([7j). The first term in 
the Model 1 potential can be written in terms of a superpotential, 

v(Q = w 2 (() + w'(0, 

W(Q = ij^j. (40) 

The first term in the Model 1 potential therefore can be viewed as one part of a super- 
symmetric partner potential system, with the second term providing explicit supersymme- 
try breaking. Since the superpotential is non-singular, the energy spectrum is necessarily 
positive-semidefinite [13] in the limit of vanishing A. 

For Q ^> 1 the potential takes the form of the potential of conformal quantum mechanics, 

V(C) = a± (41) 

with a = —1/4 — A. Therefore Model 1 can also be viewed as conformal quantum mechanics 
with regulated short distance singularity [Ml [15]. F° r a < —1/4 such models have a tower of 
an infinite number of bound states. In the limit of small A, it has been numerically verified 
for low n that the bound state spectrum takes on the characteristic geometrical form 

nir 

E n = -de vx, n G Z . (42) 

The constant C\ depends on the details of the short distance regularization. For Model 1 
this parameter is numerically determined to take the value C% = 5.02. The left panel of 
Fig. ^ shows the ground state and first excited state energy levels as a function of A. The 
right panel illustrates the emergence of a geometrical pattern in the energy eigenvalues. 

Once a large distance cut-off £a is introduced the conformal invariance is further broken. 
The energy of the negative energy states for positive values of A increases as the cut-off is 
lowered. The energy of the n th level turns positive approximately when the magnitude of its 
energy in the absence of a cut-off equals the ground state energy of an infinite square well 
of width twice (a and depth zero, that is, when 

9 

nir 7f 

C ie -sx « — (43) 
In particular, the ground state energy turns positive for a cut-off value of approximately 

Ca ~ • ( 44 ) 
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FIG. 6: Left panel: Model 1 ground state (n = 1, solid blue curve) and first excited state (n = 2, 
dashed red curve) energy E n as a function of the potential parameter A. Right panel: scaled ground 
state (n=l, solid blue curve) and first excited state (n=2, dashed red curve) energy E n = E n /e vOC 
as a function of A. The convergence of the two curves at small values of A illustrates the emergence 
of a geometrical pattern in the energy spectrum in this limit. 




FIG. 7: Left panel: Model 1 potential as a function of £ for cut-off £a = 10 and potential parameter 
A = 0.1 (solid blue curve) and corresponding ground state eigenfuction (n=l, dashed red curve). 
Right panel: ground state energy as a function of the cut-off Ca for A = 0.1. Also shown are 
the asymptotic ground state energy value in the absence of a cut-off (purple dotted line) and the 
ground state energy of an infinite square well with width equal to two times £a and depth equal 
to zero (red dashed curve). 

Below this cut-off value no negative energy states persist in the spectrum. The ground state 
wave function is displayed in the left panel of Fig. ^ for (a = 10 and A = 0.1. The right 
panel shows the dependence of the ground state energy on (\ for A = 0.1. 

Model 2: The normal mode equation describing the Atp fluctuations is equivalent to the 
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Schrodinger equation of a quantum mechanical model with potential 

1 C + 6 1 

K «> = -i(?TiF-VTT< (45) 

where in the body of the paper the supersymmetry breaking potential is zero, A = 0. The 
first term can be written in terms of a superpotential as 



v(0 = w 2 (0 + w'((), 
l l c 
C ~~ 2iTc 2 



^(0 = 7-07^- ( 4 6) 



Since the superpotential is singular, supersymmetry does not guarantee p2] that the spec- 
trum is positive-semidefinite when A vanishes. In fact, the spectrum contains a low lying, 
narrow ground state that continues to exist when A = 0. For ( » 1 the potential takes 



the conformal quantum mechanics form Eq.(41) with a = —1/4 — A. For positive values 



of A there exist a tower of an infinite number of negative energy states above the ground 
state. For small values of A, it has been numerically verified for the ground state and small 
n excited states that the bound state spectrum takes the form 

E = -0.56, 



E n = -C 2 e-7x, (47) 



with the constant Ci = 0.67, as depicted in Fig.([8j). The ground state wave function is 
displayed in the left panel of Fig. ^ for (a = 10 and A = 0.1. The right panel shows the 
dependence of the ground state energy on £a for A = 0.1. 

Note that the potentials of Models 1 and 2 can be continuously deformed into each other 
using the interpolating potential 

no = -\££r-^, (48) 

by varying the parameter a from —2 to 6. At some intermediate value of a the spectra 
of Models 1 and 2 undergo a rearrangement so that they match up with each other. In 
particular, while lowering a from 6 to —2 thus flowing from Model 2 to Model 1, the deep 
lying ground state of Model 2 shifts up and widens to become the lowest lying state of 
the geometric spectrum of Model 1, while the excited states in the geometric sector of the 
spectrum of Model 2 must widen to match up with the corresponding exited states in Model 
1. This is readily apparent for the flow of the ground state as seen in Fig.(|9]) as it evolves 
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FIG. 8: Left panel: Model 2 ground state (n = 0, solid blue curve), first excited state (n = 1, 
dashed red curve), and second excited state (n = 2, dotted black curve) energy E n as a function 
of the potential parameter A. Right panel: scaled first excited state (n=l, dashed red curve) 
and second excited state (n=2, dashed red curve) energy E n = E n /e ^ as a function of A. The 
convergence of the two curves at small values of A illustrates the emergence of a geometrical pattern 
in one sector of the energy spectrum in this limit. 



into the ground state as seen in Fig.(j7]). For A = 0, the groundstate energy approaches zero 
exactly as the parameter a approaches the value —2 from above according to 

E = -Ce~^, (49) 



for values of a near —2, as shown in the right panel of Fig. ( 10 ) . The constant C is numerically 
determined to have the value C = 0.743. While Model 1 (a = —2) does not exhibit a bound 
state for A = 0, a bound state does exist for a > —2. 
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